Dynamical, dielectric, and elastic properties of GeTe 
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The dynamical, dielectric and elastic properties of GeTe, a ferroelectric material in its low temper- 
ature rhomb ohedral phase, have been investigated using first-principles density functional theory. 
We report the electronic energy bands, phonon dispersion curves, electronic and low frequency di- 
electric tensors, infra-red reflectivity. Born effective charges, elastic and piezoelectric tensors and 
compare them with the existing theoretical and experimental results, as well as with similar quan- 
tities available for other ferroelectric materials, when appropriate. 

PACS numbers: 77.84.-s,77.22.-d,62.20.D-,77.65.-j 
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I. INTRODUCTION 

GeTe is an interesting material from both academic 
and industrial perspectives. When alloyed with an- 
timony, the electronic and optical properties of GeTe 
get dramatically modified, due to change of the micro- 
scopic structure from the crystalline to the amorphous 
phascjii^ii^ This makes it a crucial base material in phase 
change alloys used in optical storage rewritable CDs and 
DVDs. 

Besides this technological interest, GeTe attracts more 
academic oriented interest for its ferroelectric properties. 
At higher temperatures, it possesses the highly symmet- 
ric , paraelectric, rocksalt cubic structure (space group 
F3m No. 225). Below a critical temperature Tc, it stabi- 
lizes in a lower symmetry ferroelectric structure (space 
group R3m No. 160) with Ge and Te ions being displaced 
from ideal rocksalt sites. The ferroelectric transition is 
characterised by the softening of a zone-centre transverse 
optic (TO) phonon mode propagating in the [1 1 1] direc- 
tion, and the freezing-in of a relative displacement of the 
crystal sublattices4 Unlike other IV- VI telluride based 
materials such as SnTe and PbTe which have very low 
Tc being ~ 140 K and 2 K respectively, GeTe has Tc ~ 
720 K, which makes it the simplest ferroelectric material 
existing at room temperature, with only two atoms per 
primitive cell. 

Due to its interesting properties as ferroelectric and 
phase change material, GeTe has been the subject of 
many experimental and theoretical studies. The elec- 
tronic, structural and optical properties have been in- 
vestigated in the different crystalline, liquid or amor- 
phous phases Thc cublc phase instabil- 
ity and pressure induced phase transition has been the is- 
sue of discussion in some studies.— i^'-'^^i-'^-'^'-'^^i"'^^'"'^'^ However 
the dynamical, mechanical and piezoelectric properties 
of GeTe have been largely left aside. 

A Raman inelastic scattering study of GeTe was car- 
ried out in the early work of Steigmeier et al.^ who found 



two principal peaks, a first peak at a frequency of 98 
cm"^, attributed to a degenerate E mode, and a second 
peak at 140 cm~^, attributed to a Ai mode. A more 
recent experimental study, has reported values of 80 and 
122 cm~^ for the E and Ai modes respectivelyjS. 

Calculations of dynamical properties of GeTe, using 
density functional perturbation theory (DFPT) were re- 
ported by Zein et al.^ who studied the GeTe in the high 
temperature rocksalt structure. They found soft phonon 
at F and reported a value of ~ 10. 2e for the Born effec- 
tive charge. Dielectric properties of the rocksalt structure 
have also been investigated by Waghmare et al.^ who 
reported the Born effective charge, LO-TO splitting, and 
the optical dielectric constants. However, the dielectric 
properties of GeTe in the stable ferroelectric phase have 
not been examined. Ciucivara et al^ used ah initio cal- 
culations reporting a value of 10. lie for the Born effective 
charge in the ferroelectric phase. Whereas such a value is 
very close to the value previously reported for the rock- 
salt structure^ it is in clear contrast with the behaviour 
for perovskites where large Born effective charge mod- 
ifications due to the ferroelectric transition have been 
reported. 

The aim of our work was to perform a comprehensive 
study of the electronic, dynamical, dielectric, elastic and 
piezoelectric properties of GeTe in its low-temperature 
rhombohedral phase. Doing so, we uncover some prob- 
lems with previous calculations. The correct understand- 
ing of the bulk ferroelectric phase of GeTe is a prerequi- 
site for future investigations of the ferroelectric properties 
of the more complex GeTe nanostructuresJii^ 

This paper is organized as follows: In Sec. II, we 
present the details of the methods used in the present 
study. In Sec. HI, we discuss the ground-state structural 
and electronic properties of GeTe. In Sec. IV, we present 
the calculated Born effective charge and optical dielectric 
tensors. In Sec. V and VI, the phonon band structure 
and IR reflectivity are presented. In Sec. VII and VIII, 
we present the elastic and piezoelectric properties. 



II. TECHNICAL INFORMATIONS 

All the calculations have been performed using plane 
waves and norm-conserving pseudopotentials, as imple- 
mented in the ABINIT code^^ The dielectric, dy- 
namical, elastic and piezoelectric properties have been 
evaluated within the density-functional perturbation 
theoryi^i^i22ii^ We employed Hartwigsen-Goedecker- 
Hutter pseudopotentials generated including spin- 
orbit coupling, within the local-density approximation 
adopting the Teter Pade parameterizationi^Although 
this approximation is of frequent use, and gives correct 
trends, for the study of dielectric materials, one should 
be aware of some inherent limitations, due to the neglect 
of the polarisation dependence, present in the "exact" 
functional The inclusion of spin-orbit coupling does 
not lead to noticeable changes for most of the results 
presented here, except for the electronic band structure, 
which was already knowni^ Single-particle wave func- 
tions were expanded using a plane wave basis up to a 
kinetic energy cutoff equal to 15 Ha. The Brillouin zone 
integration was performed using special fc-points sampled 
within the Monkhorst-Pack schemci^ Even though the 
electronic and structural properties were found to be well 
converged using 4 x 4 x 4 fc point mesh, a denser mesh 
of 16 X 16 X 16 fc-points was required to describe well the 
vibrational properties. For reasons of consistency, the 
latter grid was used throughout this work for all calcula- 
tions. 
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FIG. 1: Calculated band structure of GeTe 

(L, (0.35,0.2,0.2) ) is 0.28 eV slightly bigger than the tun- 
neling spectroscopy results 0.2 cVM- Since the electronic 
gap is quite sensitive to the distortions, it was suggested 
by Rabc et al. that such overestimation of the gap can be 
due to the difference between the theoretical structural 
parameters used in the calculations and those of the thin 
films used in the tunneling spectroscopy measurements.— 

IV. DIELECTRIC PROPERTIES 



III. STRUCTURAL AND ELECTRONIC 
PROPERTIES 

The ground state rhombohedral structure of GeTe 
(space group R3m) has been represented in our cal- 
culations taking the z axis along the [111] of the 
conventional distorted rocksalt structure, with prim- 
itive translation vectors (f;~j^;f)i (*^' iTa' t-* ^^'^ 
(— |, —j^, f ), where the present a and c are related to 
ao and a of Table U via the relations a ~ 2aosin(^) and 
c = flo \/3(l -I- 2cos(a)). In Table U we report the calcu- 
lated structural parameters. The lattice parameters ao, 
the angle a, and the shift of the Ge and Te sublattices 
arc reported using the distorted rocksalt setting : atomic 
positions are Ge (0 0) and (0.5-t 0.5-t 0.5-r) for Ge 
and Te respectively. Our results are globally in good 
agreement with previously reported ah initio results.—. 
Moreover the deviation between the calculated structural 
parameters(ao,Q;) and experiment is less than 2%. 

In Figure [1] we show the energy band structure cal- 
culated using the theoretical structural parameters from 
Table [TJ The general features of the band structure agree 
with previously reported ah initio calculations.— The 
calculated direct energy gap at L is 0.48 eV, this value is 
slightly larger than the previously reported ah initio re- 
sults of 0.401 and 0.369 ^ The indirect gap exists at 



We have first calculated the Born effective charge ten- 
sor of atom K, defined as the induced polarization of the 
solid along the direction i by a unit displacement in the 
direction j of the sublattice of atom k at vanishing elec- 
tric field E 

dP 

Using Eq. ([1]), the Born effective charge tensor can 



TABLE I: Calculated structural parameters of GeTe. The 
lattice parameter ao (in A), the angle a (in deg), the deviation 
from of the Ge sublattice from the 0.5 sublattice position r, 
and the volume Q, = (ao/4)sin^a (in A^). The calculated 
lattice constants a and c (in A) of the equivalent hexagonal 
representation of the unit cell are also presented. 



ao a r f2 a c 



Present 5.893 88.96 0.0236 51.15 8.258 10.391 

Theoryi^ 5.886 89.24 0.0217 50.96 

Expi^ 5.98 88.35 0.0248 53.31 

Exp^Sa 5.996 88.18 0.026 
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be calculated using DFPT or within the finite-difference 
method (FDM) . In the latter case the polarization is cal- 
culated using the Berry phase technique.— We have em- 
ployed both methods in this study and the results are 
shown in Table [TTl Due to the symmetry properties of 
RSm structure, the Born effective charge tensors of Ge 
and Te are diagonal, with two independent components, 
along (Z|p and perpendicular (Z^) to the trigonal axis. 



TABLE II: Components of the Born effective charge tensor 
as calculated by DFPT and FDM, and the calculated optical 
dielectric tensor in units of vacuum permittivity eo- 



DFPT 


FDM 


zi z; 


^± ^11 £± £ll 


6.89 4.69 


6.99 4.72 59.43 49.23 



The Born effective charge tensor is strongly anisotropic 
with a difference of 2.2 between Z*^ and Z'^ components. 
Both Z'^ and Z*^ are significantly larger than the nom- 
inal ionic value of -1-2 for Ge and —2 of Tci^ This is a 
consequence of partial hybridization of the p orbitals of 
both Ge and Te which results in a mixed ionic-covalent 
bond. 

Our results for the Born effective charge of the fer- 
roelectric phase differs widely from the previous results 
obtained using the Berry phase and FDMii^ A value of 
= 10.11 was reported which is much larger than our 
calculated value of 4.72. However, the agreement be- 
tween Z* from DFPT and FDM as clearly seen in Ta- 
ble [Til q uestions the accuracy of the results obtained in 

The values of both Zy and Z'\_ are also smaller than 
those theoretically obtained for rocksalt high tempera- 
ture phasci^ii^ It has been noted previously that, for 
perovskite systems, Z* is strongly dependent on the ge- 
ometryJ^ To get more insight about the evolution of Z* 
due to the paraelectric-ferroelectric phase transition, we 
have calculated Z* in the rocksalt structure, then man- 
ually displaced atoms toward the ferroelectric positions 
(A) with small steps (St) calculating Z* at each step, 
keeping the unit cell structure fixed at the high temper- 
ature F3m lattice vectors. 

As seen in Figure [D , the Born effective charge in the 
rocksalt phase is 10.68, in good agreement with the thcor- 
tical results of Waghmare et al.^ who reported a value 
of 10.8. It has been shown previously that such high 
value of Born effective charge can be due to extremely 
narrow electronic gap possessed by the high temperature 
rocksalt phased We have noticed that both Z]_ and Zj* 
drop as the atoms approach the ferroelectric sites. Such 
modification of charge with atomic position is considered 
a feature similar to what was obtained previously in case 
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FIG. 2: (color online) Evolution of longitudinal Born effec- 
tive charge Z*^ (blue) and transverse Born effective charge 
(red) as a function of atomic displacement from the rocksalt 
ideal positions. 

of ABO3 perovskites.— 

Similar to the effective charge tensor, the calculated 
dielectric tensors arc diagonal consisting of two indepen- 
dent components, parallel and perpendicular to the trigo- 
nal axis. Our value is larger than the previously reported 
experimental value (35-37.5).— Generally speaking the 
high value of (£°°) comes as a consequence of the low 
value of the electronic gap. 

V. DYNAMICAL PROPERTIES 

Since there are two atoms per primitive unit cell, there 
will be six phonon branches. The phonon branches are 
divided into three acoustic and three optical phonon 
modes. Along the T — T direction (trigonal axis), the 
phonon dispersion curves can be classified as E or Ai, 
according to whether the atomic displacements are per- 




FIG. 3: Calculated phonon dispersion curves of ferroelectric 
GeTe at the theoretical lattice parameters. 
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pcndicular or parallel to the trigonal axis, respectively. 
Along the other directions, the branches cannot be clas- 
sified as pure E or Ai modes. 

In Figure[31 we show the calculated phonon band struc- 
ture obtained by DFPT. The general feature of the band 
structure are close to those of Bi which has a similar 
structural unit cell.— Because of the nonvanishing com- 
ponents of the Born effective charge tensors, the dipole- 
dipole interaction had to be properly included in the cal- 
culation of the interatomic force constants.— Such inclu- 
sion of the dipolc-dipole interaction in the interatomic 
force constants results in splitting of longitudinal and 
transverse optic modes (LO-TO splitting). Within this 
treatment, the Ai and E modes (for wavevector aligned 
with the trigonal direction) are characterized as LO and 
TO modes respectively. The crystal being uniaxial, the 
LO-TO splitting at T might vary with the limiting direc- 
tion. We actually notice a very strong dependence of the 
LO-TO splitting on the wavevector. 

It is difficult to grow large high quality GeTe crys- 
tals, so, little experimental information on its vibrational 
modes has been collected. The comparison between the 
calculated mode frequencies at T with those measured 
by Raman scattering experiments is also difficult j^i^ It is 
quite known that GeTe samples are non-stoechiometric 
and contain a high concentration of free holes.— This re- 
sults in complete screening of the dipole-dipole interac- 
tion by the conduction electrons at zone centeri^ How- 
ever, the comparison can be possible if we artificially en- 
force such a complete screening of the dipole-dipole in- 
teraction at the zone center. The results are shown in 
Table IIIIl the first row contain the values of frequencies 
as calculated by including the dipole-dipole interaction 
using a semiconducting screening, the second row con- 
tained the values of frequency as calculated by complete 
screening of the dipole-dipole interaction.— We note that 
the values calculated by complete screening of the long 
range force are in more agreement with the existing Ra- 
man inelastic scattering resultsi^i^ Note that in absence 
of the LO-TO splitting as a result of complete screen- 
ing of dipolc-dipole interaction, the Ai mode has been 
characterized experimentally as a TO mode. In fact, 
it was suggested that the softening of this mode at the 
critical temperature is responsible for the ferroelcctric- 
paraelectric transition.— 122 



VI. IR SPECTROSCOPY 

The dielectric tensor e{uj) in the lowest frequency range 
can be related experimentally to the IR spectra. It can 
be calculated theoretically by accounting for ionic re- 
laxations in the calculations of the permittivity tensor. 
The ionic contribution to e(uj) comes mainly from opti- 
cal phonon contributions (without damping) to the IR 
oscillator strength Sm,ij 



TABLE III: Phonon frequencies at the zone center (in cm ^) 
calculated for a q vector || to trigonal axis. 



EjTO) AijLO) 



Semiconducting 71 153 

Complete screening 100 134 

Exp^ 98 140 

Exp^ 80 122 



. TO ^ 

m 

where is the volume of unit cell and m is the phonon 
mode rank. 

In Table IIVI we present the calculated IR oscillator 
strength and the components of e(0) for a perfect semi- 
conductor GeTe crystal. From Table IIVI it is obvious 
that the ionic contribution to £{0)± comes purely from 
E modes and to £(0)|| from Ai mode. The inclusion of 
the ionic contribution results in a strongly anisotropic 
e(0) with e(0)_L almost three times larger than e(0)||. In 
fact, the strong anisotropy of e(0) has been expected in 
regard of the frequency difference between E{71 cm^^) 
and Ai(119 cm~^) modes. Another reason for the large 
anisotropy comes from the fact that one of S± of E modes 
is almost twice S^^ of Ai mode. 

In Fig. [3]wc present the calculated IR reflectivity as- 
sociated with E and Ai modes. The reflectivity related 
to Ai mode can be associated with light incident parallel 
to the trigonal axis of the crystal i.e., perpendicular to 
(001) surface. Similarly, the E modes reflectivity shall 
be associated with light incident perpendicular to the 
trigonal axis, i.e., parallel to a (6) directions which are 
perpendicular to 100 (010) surfaces. 

The range of light wavelength at which the maximum 
reflectivity occurs differs depending on the surface. For 
example, whereas the maximum reflectivity for light inci- 
dent on (GDI) surface occurs for wavelengths in the range 



TABLE IV: Calculated values of IR oscillator strength tensor 
5 ( X 10"'') (in a.u.) of optical modes and components of static 
permittivity tensor in units of eo- 



i II 


S{E) 


5.74, 50.45 


0.00 


S{Ai) 


0.00 


26.04 


e{0) 


254.89 


81.58 
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FIG. 4: (Color online) Calculated IR reflectivity spectra of 
GeTe (without anharmonic damping). 



65-83 /im, it is in the range of 68-137 /xm in case of 100 
or 010 surfaces. 



VII. ELASTIC PROPERTIES 

In this section we present our results for the elastic and 
compliance tensors. The elastic tensor is defined as the 
change in the stress of the solid in the direction a as the 
strain is changing in the direction /3, where a,/? = 1...6 
in Voigt notation 



Cq/3 



drip 



(3) 



The above equation split into two main contributions 



dag 
drip 



E 



darr du 



k.i 



duk/ drip 



(4) 



The first term is the frozen (clamped) ion elastic tensor 
(c°), the second term includes contributions from force- 
response internal stress and displacement-response inter- 
nal strain tensors. The second term accounts for the ionic 
relaxations in response to strain perturbations. The ad- 
dition of the two contributions is the relaxed ion elastic 
tensor c. The compliance tensors is simply defined as the 
inverse of the elastic tensor. 

In Table El we show our results for the elastic and 
compliance tensors. Due to the low symmetry of R3m 
of rhombohedral phase of GeTe, there are six indepen- 
dent elastic constants. The obtained values for the elas- 
tic tensor constants satisfy the mechanical stability re- 
strictions for trigonal type unit cells cn — \ci2\ > 0, 

(Cll -I- C12)C13 - 2c^3 > 0, (Cii - C12)C44 - 20^4 > 04^ 

Usually, the inclusion of internal relaxations reduces the 
values of elastic tensor components due to the relief of 
stress tensor components. The same behavior is also 
reflected in the increase of the compliance tensor com- 
ponents. The effect of internal relaxation is more pro- 
nounced in the case of C33 which represents the axial 



shear along the surface perpendicular to the three-fold 
rotation axis. However, one should emphasize here that 
the effect of internal atomic relaxation is not as strong as 
for other materials such as ZnO and BaTiOai^ The small 
differences between the clamped and relaxed compliance 
tensor might explain the similarity of the bulk modulus 
of the ferroelectric phase and the rocksalt phase as has 
been suggested by the experimental findingsJ^ 



TABLE V: Clamped ion c° and relaxed ion c elastic tensors 
(GPa), clamped ion s° and relaxed ion s compliance tensors 
(in TPa-M. 



index 




c 




s 


11 


116.16 


112.74 


11.66 


12.28 


12 


22.34 


19.84 


-2.35 


-2.92 


13 


36.78 


27.56 


-3.94 


-4.31 


14 


27.08 


24.68 


-5.82 


-8.45 


33 


86.98 


59.88 


14.83 


20.66 


44 


65.29 


44.37 


20.14 


31.93 



The bulk modulus can be readily calculated, using the 
above results, from the compliance tensor )^ 



B 



3 

Q/3 



(5) 



Using the above equation, we get a value oi B = 45.17 
GPa which is in quite good agreement with the experi- 
mental value of 49.96±3.2 GPaJ^ 

We have also calculated the bulk modulus by fitting 
the total energy as a function of volume curve using the 
Murnaghan equation of state. The obtained value of i? = 
44.3 GPa is in much better agreement with experimental 
findings than the theoretical value reported in Ref. ITsj . 

Unfortunately we are not aware of any experimental 
result for the elastic constants of GeTe. The agreement 
between the value of the bulk modulus value extracted 
from the elastic tensor and those calculated by fitting 
to the equation of state or measured in experiment nev- 
ertheless insures the overall reliability of the calculated 
values of the elastic tensor. 



VIII. PIEZOELECTRIC PROPERTIES 

The proper pizeoelectric tensor e is defined as the in- 
duced polarization in i direction due to a strain change 
for index a 



dr]a 



(6) 



6 



dP' 



dr]a 



E 



dP/ du. 



•ik 



(7) 



Where P'i is the reduced (rescaled) polarization as de- 
fined in the appendix of Ref. H^- The first term is the 
proper homogeneous strain contribution to the piezoelec- 
tric tensor e" which arises mainly from the sole elec- 
tronic contribution. The second term, often called in- 
ternal strain piezoelectric tensor, includes contributions 
from the Born effective charge tensor diagonal and in- 
ternal relaxation. The second term represents the ionic 
contribution to piezoelectric tensor. 



TABLE VI: Independent components of proper homogeneous 
piezoelectric tensor e", internal strain piezoelectric tensor, 
and the proper total piezoelectric tensor (in C/m^). The 
piezoelectric constant tensor d (in pC/N) and the electrome- 
chanical coupling constants k are also given. 



Similar to the zero strain dielectric tensor e, the stress 
free dielectric tensor e"^ is diagonal, with two indepen- 
dent values £j[ ^ 90.70e" and e3;=360.09e". Note e" have 
components larger than s as expected for any piezoelec- 
tric material. 

In the last two rows of Table IVIl wc show our cal- 
culated piezoelectric constants d and electromechanical 
coupling constants. Even though the various components 
of d are quite smaller than those reported for PMNT sys- 
tem;^ it is interesting that our calculated value of c?i5 is 
comparable with that of the giant piezoelectric materials 
PMN-PT and PZN-PTi^ where dis extends between 131- 
190 pC/N. On the other hand, the calculated electrome- 
chanical coupling factors in general are far less compared 
with that of PMN-PT and PZN-PTi^ or PMNTi^, how- 
ever they are slightly better than those of ZnO^ 





15 


21 


31 


33 




0.88 


1.24 


0.89 


-0.31 


internal strain 


-6.01 


-0.69 


-0.95 


-2.78 


Total 


-5.13 


0.55 


-0.06 


-3.09 


d 


-173.00 


51.68 


12.78 


-63.44 


k 


0.54 


0.26 


0.12 


0.49 



In Table IVIl we present the results for the piezo- 
electric tensor. There are 4 independent components. 
The components 631 , 633 represent the induced polariza- 
tion along the trigonal axis created in response to shear 
strain in the a6-plane and along the trigonal axis respec- 
tively. The other components describe the induced po- 
larization along the primitive axes (a, 6) by shear strain. 
The calculated value of 633 is -3.09 (C/m^). The proper 
homogeneous strain contribution has been found to be 
-0.31(C/m^), however the strain contribution is much 
larger, adding -2.78 (C/m^). The strain contribution was 
relatively large in the case of 615 with a value of -6.01 
(C/m^). The strain contribution of the last case being 
twice the value of 633 can be explained by the noticeable 
anisotropy in the Born effective charge, and by the large 
strained-induced ionic motion in the lateral direction in 
response of a strain applied along the trigonal axis. For 
the other elements, the strain contribution reduces the 
polarization, which is almost cancelled, as in the case of 
631- 

The efficiency of the produced electric energy vs the 
spent mechanical energy can be estimated by calculating 
the electromechanical coupling constant, defined as 



where dia is the piezoelectric constant dia = Sapeip, and 
e'^ is the free stress dielectric tensor, related to the vanish- 
ing strain dielectric tensor shown in Table IVIl via Eq. 20 
of Ref. m . 



IX. CONCLUSION 

We have investigated the dielectric, dynamical and me- 
chanical properties of the ferroelectric phase of GeTe 
within density-functional perturbation theory. Our study 
covers all the linear couplings between applied static ho- 
mogeneous electric field, strain, and periodic atomic dis- 
placements : Born effective charge, dynamical matrix at 
the zone-center, clamped and dressed elastic constant, 
optical dielectric tensor, adiabatic dielectric tensor, free 
stress dielectric tensor, piezoelectric coefficients, elastic 
and compliance tensor. We also examined the phonon 
band structure, as well as the change of the Born effec- 
tive charge tensor with the atomic positions. Concerning 
the latter, our results disagree with those of Ciucivara 
et alt^ by a factor of two. Since we observe a rapid 
variation of the Born effective charge with the atomic 
positions, we hypothetize that their result has not been 
obtained for the ferroelectric phase, but for the symmet- 
ric paraelcctric phase. Our other results have also been 
discussed, and compared with the available theoretical 
and experimental results. 
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